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ABSTRACT 


This report summarizes the progress achieved to data under NASA Grant 
NSG 1213. In Section I, some observations on the perfect noninteraction 
for a discrete parameter set is given. Section II investigates the sensi- 
tivity of a decoupled flight control system with respect to system parameter 
variations. In Section III, a brief description of a set of computer 
programs developed for this project is given. The program listings will 
be furnished upon r^uest. 


I. Perfect noninteraction for a discrete parameter set 


When parameter variations are defined over a discrete set of parameter 
values, it may be possible to compute a fixed control law which yields 
decoupled behavior for all admissible parameter values. This possibility 
was pursued as a first step to the decoupling problem since when a solution 
exists, it may be obtained by some fairly simple ccanputations. 

It has been shown by Gilbert [l] that the class of all feedback control 
laws of the form, 

u = Fx + Gv, 

which decouple the system, 

i = Ax + Bu, y = Cx, 

must have the following structure for F and G, 
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where the matrices A ,D,J, and G, and the scalers TT„ are determined from 

the given matrices A,B and C, while and are free design parameters. 

In addition a determine the poles of the closed-loop decoupled system, 

Jk 

while X . determine the scalar gain factor of each diagonal element in the 

closed- loop transfer matrix. For each system parameter set A., B. , there 

i i ^ ^ 

is a corresponding F., G. with free parameter and X . The idea is then 

^3- X ^ Jk j 

to seek values of a., and X . such that 

Jk 0 


Fi = Fa = F3 — 


Gi = Ga = Gs = 

This yields a system of linear equations in the free parameters. If a solu- 
tion to these linear equations exists, then a fixed control law exists 
which decouples the system for all possible values of A, B and C. Con^juter 
programs are being developed for the generation of the appropriate linear 
equations and their solution. If no solution exists, perfect decoupling 
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is not possible over the given parameter set and one must proceed to a more 
complex design algorithm, such as the guaranteed-cost procedure described 
in the original prt^osal. The requirement of stability of the closed-loop 
system places further constraints on the existence of a solution. Finally 
the requirement of invariant closed-loop poles may be imposed by requii^jig 
that the same values of and be used in each and G^. However 
the places even more constraints on the existence of solutions. 

An abstract geometric approach to the problem of perfect decoupling 
over a discrete parameter set has recently been developed by Asher and 
Mulholland in C2]. However the approach presented here is more direct arwi 
easier to compute. 
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II. Sensitivity of decoupled system 

The following example is taken from Pope [3]* and involves the lineari- 
zed longitudinal equations for a landing maneuver ( STOL aircraft with a 
trim speed of 60 knots and a 7^ degree angle of attack). Let x denote the 
syston state with cos^onents, 
xi = forward speed 
X2 = pitch angle 
X3 = pitch rate 
X4 = vertical velocity 

and let u denote the control vector with components, 
ui = elevator deflection 
via = thrust change. 

Then the system dynamics are given by (see page 23 of [3]) 


X = Ax + Bu, y = Cx, 


( 1 ) 


where 


A = 


-.032 -32.2 0 .133 

0 0 10 
.00137 0 -,743 -.0014 

-.02 4.2 96.5 -.3 


( 2 ) 


0 .65E-3 


B = 


0 0 

-.989 -.7E-5 



0 1 
0 0 


0 

0 


(3) 


3.0 -.87E-3 
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In [2] a state feedback control law of the form u = Fx + Gv is 

determined which decouples pitch rate and forward speed. The decoupling 

values of matrices F and Q are computed from computer programs developed 

by Gilbert and Pivichny [U], For the nominal parameter values given in 

(2) and (3)> the following values are computed for F and G, see page 40 of 
- -* 

C3J , 


.01193 .6605 .8665 .3E-4 


F = 


-1.4892E3 4.9538E4 0 


-204.615 


G = 


-.088 -.0229 

0 3.2308E3 


(4) 


(5) 


A computer program was developed for the computation of the closed- 
loop transfer matrix H(s), with elements h^^(s). The computer program was 
used to verify that the matrices given by (4) and (5) did indeed decouple 
the nominal system. The theory was verified with diagonal elements, 




0.087s 

1^ + 1.6s + 1 


hggCs) 


2.1 

s + 1 


as predicted in [3]. Off diagonal elements, while not identically zero, 
were smaller than the diagonal elements by a factor of IE-5. The closed- 
loop characteristic polynomial is (s® + 1.6s + l)(s + l)(s + .12199) • 

However contrary to the claims in [3], the solution appears to be very 
sensitive to parameter variations. For example, if the system parameter 
values are rounded off to two significant figures and the same values of 
F and G are used, the tremsfer element Ifei(s), which couples elevator 
deflection to forward speed, is no longer negligibly small. In particular, 
for the parameter variations corresponding to the above round-off modifica- 
tion of system parameters. 


* The values of F and G given on page 40 are evidently in error. Correct 
values are given here. 
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b>i (s) 


-.039?(s® + 70.15s + 19.75) 
s* + 1.761s® + 1.098s® - .0973s - .I65U 


s 


The coupling between thrust and pitch rate remained essentially zero (iha(s) 
was smaller than raaaining elements by a factor of lE-5)* The closed-loop 
system is unstable, with a pole at s = 0 . 33 * 

If the system parameters are modified to reflect a new trim speed of 
55 knots, the A and B matrices become* (page 88 of C3])» 


A = 


• .02912 

-32.2 

0 

.12103 

) 

0 

1 

0 

00125 

0 

-,676 

-.00127 

-.0182 

4.2 

89.829 

-.273 


0 .55E-3 


B = 


0 0 

-.831 -. 6 e -5 


L 2.52 -.73E-3J 

This corresponds to parameter variations of approximately 10 ^, If the 
nominal values of F and G are used for this new flight condition, decoupling 
between elevator deflection and forward speed is again destroyed. In this 
case hai(s) beccanes. 


_ -.0397(8® 70.15s + 19.75) 

s* + 1.761s® + 1.098s® - 0.0973s - 0.1654 


* There is apparently a sign error in the elements a*i and b42 listed on 
page 88 . 
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As in the previous case^ parameter variations result In aii unstable closed* 
loop system^ In this case the unstable closed* loop pole Is located at 
s s 0.28. This numerical example will be used to test the design procedure 
developed In this study. 
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III. Ccagputer programs 


In this section, we describe the set of computer programs developed 
for this research pro;]ect. This set of c<miputer programs falls into two 
categories (i) a digital simulation package and (ii) numerxcal optimization 
programs. All these programs are written in standard Fortran IV language. 
Hence they can be run cai almost any computer installations. Some of these 
programs could be very time consumming when applied to hlgh>order systems. 

It is planned to run most of them on Hewlett-Packard 2100 minicomputer 
which is owned by the lAiiversity of Colorado at Colorado rings. With 
proper arrangement of the (iterating system of the mlnicompi. ter, these 
program can be running in an interactive mode. 

A brief description of these progrms is given below. 

(l) Digital simulation package 

This program is baaed on the Runge-Kutta fourth-order algorithm. 

The input data to the program is the parameter of the plant, observer, noises 
feedback and feedforwen^i control laws, command input v, initial states 
Xo , , initial and final time to and ti . The over-all system configuration 

is shown below. The output of this program is the time history of the system 
outputs, which can be either printed in the graph form or stored in an array 
to be used in conjunction with the optimization programs described later. 
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(ii) Optimization programs 

Two programs are Inclwled here, Box Complex Algorithm and Rosenbrock 
Hill Algorithm. Both algorithms find the maximum or minimum of a multi- 
variable, nonlinear function subject to nonlinear inequality constraints: 
Optimize F(X^, Xg, ... , X^j) 

Subject to \ \ ^ \ » k = 1, 2, ..., M 


The implicit variables •••> dependent functions of the 

explicit independent variables, X^,Xg, Xj^. The upper and lower 

constraints and are either constants or functions of the independent 
veurlables , 


The Box Complex Algorithm is a sequential search technique which has 
proven effective in solving problems with nonlinear objective functions 
subject to nonlinear inequality constraints. No derivatives are required. 
The procedure should tend to find the global maximum due to the fact that 
the initial set of points are randomly scattered throughout the feasible 
region. The Rosenbrock Hill Algorithm proceeds per the imcons trained 
Rosenbrock procedure until convergence is reached or a boundary zone in the 
vicinity of the constraints is entered. Both of these two programs are 
taken from [ 5 ], which contains detailed Information. 
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